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ABSTRACT 

The numerical study of QED radiative corrections in the low x region for the 
HI experiment at HERA was performed. The analytical program TERAD91 
was used to obtain the values of the radiative correction factor for the first 
measurement of the proton structure function F 2 (x, Q 2 ) in the range x = 10~ 2 — 
10" 4 and 5 < Q 2 < 80 GeV 2 . 

In order to extract from experimental data the inclusive cross section of deep 
inelastic ep-scattering 

e(k 1 )+p(p 1 )^e{k 2 )+X{p 2 ) (1) 

in one-boson-exchange approximation one must to take into account contributions 
from higher order QED 1 and electroweak processes. For example the bremsstrahlung 
process: 

e(h) + p( Pl ) - e(k 2 ) + X(p 2 ) + 7 (fc), (2) 
and also the elastic radiative tail: 

e(fci) + p(pi) -> e(k 2 ) + p(p 2 ) + 7 (fc). (3) 

contribute to the observed cross section of the reaction (1). 

In the long period from the SLAC experiments 2 till the operation of HERA, the 
experiments with neutral current deep inelastic ep-scattering had a relatively simple 
structure, and the measurements completely rested on the registration of energy and 
angle of the scattered electron. The new detector generation, the HERA detectors 
HI 3 and ZEUS, 4 allows to measure both the scattered electron and hadronic final 
state. 

For the physical analysis of the ep-collisions at HERA one can use not only the 
familiar electron variables 



Ql = -{h- h) 2 , y e =Pi(k 1 -k 2 )/p 1 k 1 , x e = Ql/(sy e ), (4) 



where 

s = (h + Pl ) 2 = AE e E p , (5) 

but also the kinematical variables from the hadron measurement 

Ql = -(P2 -Pi) 2 , yh=Pi(p2-pi)/pih, x h = Q 2 h /(sy h ), (6) 

or some mixture of both. Here E e and E p are the energies of incident electron and 
proton. 

For different choices of variables there are differences in the predictions for the 
contribution of the higher order processes. 

The theoretical investigations of these radiative effects, or radiative corrections 
(RC), for experiments at HERA were summarized in Proceedings of the Workshop 
"Physics at HERA" . 5 Different programs to calculate the radiative corrections have 
been cross-checked . 6 

The analysis of neutral current deep inelastic scattering data in the HI experi- 
ment at HERA in 1992 7 was based on two independent approaches: 

i) a " detector oriented" approach ( analysis I ) 

ii) a "physics oriented" approach ( analysis II ) 

Within the analysis I the kinematics of the inclusive deep inelastic scattering 
process (1) were determined from the energy E' e of the scattered electron, and its 
polar angle 9 e , measured relative to the proton beam direction 

Ql = 4E e E' e cos 2 (6 e /2), (7) 
y e = l-(E' e /E e )sm\9 e /2). (8) 

The cross sections, acceptances of detector, radiative corrections and resolution 
effects were obtained in (^J~E~l,6 e ) bins. The cross sections, that have been measured 
using this binning, were transformed into cross sections in (x e , Qe)-bins. 

In the analysis II the kinematics of the process (1) were obtained from a 
measurement of the electron variables E' e , 9 e and scaling variable yh from the hadrons. 
The hadronic variable yh (6) can be determined using the Jacquet-Blondel method 8 

EEh ~ Pz,h / n \ 
~^2E ' ^' 

hadrons e 

where Eh is the energy of a hadron and p Zj h its momentum component along the 
incident proton direction. 

The momentum transfer Q 2 was determined from the electron variables 9 e and 
E' e as in (7) and Bjorken x by using a mixture of the electron and hadron measure- 
ments 9 

x m = Q 2 J{sy h ). (10) 



In the analysis II cross sections, radiative corrections and efficiencies were cal- 
culated in (x m , Ql)-bins. 

A further difference between the two methods of the analysis of the data was 
the use of two different types of programs — the Monte Carlo event generators HER- 
ACLES 10 in the analysis I and the analytical program TERAD91 11 in the analysis 
II — for the calculations of the radiative correction factor S(x,Q 2 ). The radiative 
correction factor 5(x, Q 2 ) is defined by the equation 

) d 2 a Bom /dxd Q2 W 

where d 2 a thcor /dxdQ 2 is the theoretical approximation to the measured cross section 
(Pa meas /dxdQ 2 , which contains the contributions from higher order electroweak pro- 
cesses, and d 2 a Born /dxdQ 2 is the Born cross section of the process (1). 

In this note we describe the numerical calculations of the RC factor 5(x,Q 2 ) 
for the "Zeuthen analysis" of the proton structure function F 2 (x, Q 2 ) in the range 
x = W- 2 - 1(T 4 and 5 < Q 2 < 80 GeV 2 in the HI experiment at HERA 7 

The Monte Carlo HERACLES is suited for the calculation of RC to the differ- 
ential cross sections in terms of the electron variables because the program calculates 
cross sections and generates events in terms of the variables x e and Q 2 . Therefore 
most of the generated radiative events in HERACLES have low Q\ and the calcula- 
tion of the radiative corrections in the terms of other variables ( hadron, mixed ) can 
be performed with less efficiency only. For this reason in the analysis II was used the 
analytical program TERAD91. 

TERAD91 - a program package to calculate radiative corrections for deep in- 
elastic ep-scattering - was created in 1991 during Workshop "Physics at HERA" and 
accumulates about 15 years experience in the field of calculation of the radiative cor- 
rections for deep inelastic /iV-scattering , 14 The program TERAD91 can be applied to 
calculate the RC factor of the order 0(a) to the measured differential cross section, 
d 2 a meas /dxdQ 2 , after the acceptance and smearing corrections. 

The program package TERAD91 contains one part, TERAD, which is based 
on the model-independent treatment 15 of the leptonic QED corrections for neutral 
current ep-scattering, including bremsstrahlung process ( 2) and QED vertex cor- 
rections. TERAD does not necessarily rely on the quark parton model and uses the 
phenomenological structure functions Fi,F 2 and F 3 , which are defined in the Born 
approximation. 

The program TERAD91 can be used for the model-independent calculations 
of QED corrections for neutral current deep inelastic ej9-scattering in the terms of 
several variables - electron, hadron, mixed variables 16 - and in the new version 17 also 
in terms of the Jacquet-Blondel variables . 18 

In the low Q 2 region the Born cross section of the deep inelastic scattering (1) 



is determined by two structure function F 2 and 2xFi = F 2 /(l + R): 



d 2 a 2ita 2 ( , , y 2 \ . , 

«? = w( 2(1 - v) + TTflJ F,(l,<n - (12) 

In order to calculate the RC one has to use a realistic parametrization of the 
structure functions Fi(x,Q 2 ) over the full range of x and Q 2 . We have used several 
different parameterizations of the deep inelastic proton structure functions , including 
MRS D- and MRS DO parameterizations . 13 For region Q 2 < 5 GeV 2 the parameter- 
izations were taken at Q 2 = 5 GeV 2 and multiplied by factor 19 (1 — exp(— aQ 2 )) , 
a = 3.37 GeV~ 2 . In our calculations we have assumed R — 0. 

In the " Zeuthen analysis" the following x, Q 2 -binning was used. There were four 
Q 2 -bins: 

Q 2 = 5 - 10, 10 - 20, 20 - 40, 40 - 80 GeV 2 (13) 
and four bins per x-decade: 

log x = -4.0, -3.75, -3.50, -3.25, -1.0. (14) 



The kinematical region of the experimental measurement was defined by: 

157.5° <6 e < 172.5°, 0.05 < y e < 0.60 (15) 
for electron variables and 



157.5 <9 e < 172.5°, y h < 0.50, 9 jet > 10° (16) 
for mixed variables. 

Figs. 1,2 show the results for the radiative corrections in the terms of the mixed 
variables x m , Q 2 . The RC in the terms of the mixed variables, S(x m , Q 2 ), are not small 
and can be reach 7 % at x = 10~ 2 . 

The "Zeuthen analysis" of the deep inelastic scattering has included also the 
analysis of the experimental data in the terms of the electron variables x e , Q 2 corre- 
sponding to the binning ( 13), ( 14). Fig.3,4 show the results of the calculations of the 
RC in the bins x e , Q 2 e) that have been obtained with the program TERAD91. 

The corrections determined by the electron variables, 5(x e , Q 2 )-, are known to be 
large at low x and dominated by hard photon emission from the incoming electron. 
From the comparison of Fig. 3 and Fig. 4 it is seen that the sensitivity of the RC to 
the choice of the structure functions is not small. 

Because the collinear radiation of photons in the final state of the reaction ( 2) 
is not resolved in the electromagnetic calorimeter, one must calculate the RC in terms 
of the calorimetric variables including the total energy of the electromagnetic cluster 



E = E' + E 

J- 1 VIS ~ J-Jry. 



(17) 



In the program TERAD91 the possibility to calculate the RC in terms of calori- 
metric variables, E vis , 9 e , has not been included yet, but with the help of the leading 
log calculations 12 one can estimate the size of the influence of the calorimetric mea- 
surement on the radiative corrections. 

The calorimetric measurement cannot distinguish between a single electron and 
an electron with accompanying photons. Therefore the final state radiation has exper- 
imentally no effect on the measurement of the kinematics of deep inelastic scattering 
and can be absorbed in the definition of the final electron energy of the process ( 1). 

This leads to the following way to take into account the calorimetric measure- 
ment. We can calculate the contribution of the final state radiation in the leading 
logarithmic approximation and subtract from the total RC, that we obtain with 
the program TERAD91. The more detailed investigation of the calorimetry problem 
composes the subject of a separate publication. 
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Fig. 1. Radiative correction factor 6(x m ,Ql) in (x,Q 2 )-h'ms for MRS D- parametrization. 





Fig. 2. Radiative correction factor 5(x m ,Q^) in (x, Q 2 )-bins for MRS DO parametrization. 
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Fig. 3. Radiative correction factor 5{x e ,Q 2 e ) in (x,Q 2 )-bins for MRS D- parametrization. 
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Fig. 4. Radiative correction factor 6(x e ,Ql) in (x, Q 2 )-bins for MRS DO parametrization. 



